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SUMMARY 


In support of the U. S. Department of Transportation's 
Climatic Impact Assessment Program, the Aeronautical Research 
Associates of Princeton have undertaken an analysis of the 
hydrogen-nitrogen-oxygen reaction systems in the lower strato- 
sphere as these are perturbed initially by individual aircraft 
engine exhaust plumes. These analyses were directed primarily 
to the question, "Do any significant chemical reactions occur 
either among exhaust chemical species or between these species 
(and their products) and the environmental ozone while the 
exhaust products are confined in intact plume segments at 
relatively high concentrations?" Secondarily, the joint effects 
of diffusive mixing and chemical kinetics on these reactions 
were investigated, using the techniques of second-order closure 
diffusion/chemistry models developed by Donaldson (ref. 1) and 
Hilst (ref. 2) at A.R.A.P. Operated in a coupled mode, these 
models permit quantitative assessments of the effects of inhomo- 
geneous mixing and diffusion limitations on the chemical reaction 
rates . 

Throughout these studies attention has been focussed on the 
larger problem of the potential depletion of ozone by S3T exhaust 
materials emitted into the lower stratosphere . Although the 
analyses completed to date are handicapped by inadequate observa- 
tional data, needed both as input to the models and for validation 
of the predictions, the following results have emerged: 

• Under a variety of conditions of aircraft emissions of NO and 
HgO and of aircraft-induced and natural turbulence, the depletion 

of Oo in the early post-vortex plume remains as a critical vari- 
able in the immediate post -vortex chemistry. Within reasonable 
emission rates of NO, a significant "O^-hole" is possible. 



•Only under what appears to be 1) rather extreme concentrations 
of water vapor and/or photolysis rates, and 2) minimum emissions of 
NO, is there a significant conversion of N0 X to HNO^ during the 
first few sunlit hours after fly-by. This conclusion could be 
tempered by more exact data on H 2 0 emissions, stratospheric turbu- 
lence, and photolysis rates, but, insofar as present estimates are 
reliable, early HNOo formation cannot be depended upon for subse- 
quent buffering of the N0 X catalytic cycle for 0^ destruction. 


INTRODUCTION 


During the past two years, the Aeronautical Research Associ- 
ates of Princeton have undertaken analyses of SST-related chemistry 
and diffusion in support of the D0T/CIAP studies. This work was 
carried on under Contracts NAS1-11433, NAS1-11873 and NAS1-12475. 
This report constitutes the final report on NAS1-12475, and has 
also been prepared for publication, in a somewhat amended form, in 
fulfillment of NASI -11873. 

The present report is one small part of the analyses which the 
U. S. Department of Transportation's Climatic Impact Assessment 
Program (CIAP) and a world-wide community of scientists are assemb- 
ling. It is also a very specialized study, directed to the pro- 
cesses of chemical reactions and turbulent diffusion in individual 
SST exhaust plumes during the first few hours of their existence. 
Prior to the CIAP analyses this phase of the overall problem was 
generally ignored, although it was clear that the genesis of any 
environmental problems must be the individual aircraft's contribu- 
tion. In particular, during the first few hours or days of their 
existence these exhaust products can be expected to be confined to 
relatively intact plume segments at concentrations much in excess 
of their background levels in the lower stratosphere. (We have 
termed this phase the "highly perturbed stratosphere.") Noting 
this fact, the question arises, "Does anything significant happen 
during this period of high concentrations?" 

The next logical question, of course, is, " Can anything 
significant happen?" Several possibilities worth exploring come 
to mind. For example, given the very stable temperature structure 
of the lower stratosphere, the exhaust materials may be confined 
in these plume segments for extended periods of time, periods 
during which they are not in chemical contact with the environ- 
mental constituents with which they may react, such as ozone. If 
this is the case, the local environmental impact of the exhaust 
materials could be severe, but these volumes could also be such a 
minute fraction of the total stratosphere that such effects would 
be totally negligible. 
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As another example , chemical reactions which convert environ- 
mentally hazardous exhaust constituents into relatively inert 
materials can proceed at a much more rapid rate when high concen- 
trations obtain. A prime example, which is explored in depth here, 
is the conversion of NO to N0 2 and thence to HNO^ in the presence 
of sunlight, O3, and H 2 0. This complex process, which would 
effectively buffer the N0 X catalytic cycle for the destruction of 
Oo, is particularly effective when the concentration of H 2 0 is well 
above background stratospheric levels. However, this process also 
depends upon an adequate amount of Oq, which must be supplied by 
either local chemical production or by turbulent entrainment of 
the neighboring environmental Oq. In order to assess this possi- 
bility we must, therefore, work^with coupled chemistry and diffu- 
sion simulation systems. 

As the reader may judge for himself from what follows, these 
studies of the chemistry and diffusion ol the highly perturbed 
stratosphere suggest that there may be important chemical processes 
which go on during the early life of individual SST plumes. How- 
ever, the question as to whether or not these processes actually do 
occur cannot be answered now. Much harder data on the small-scale 
turbulence in the stratosphere and a critical testing of the chemi- 
cal chains assumed in these analyses, particularly as they actually 
occur in the stratosphere are required. Happily, analyses such as 
the present one point the way to obtaining this critical informa- 
tion. 


THE PROBLEMS POSED 


A qualitative description of the pertinent features of an SST 
engine exhaust plume during the first few hours after a single SST 
has flown by in the lower stratosphere provides the framework for 
the quantitative treatment of the diffusion and chemistry of the 
exhaust materials and environmental constituents which follows. 
First, it is important to recognize that the lower stratosphere 
(15-22 km above mean sea level) is characterized by generally 
strong hydrostatic stability. As a result, local wind shear- 
produced turbulence is practically nonexistent (ref. 2 ) , Larger— 
scale meandering wind currents, and some _ vertical wind velocity 
variations are observed, but the production of small-scale turbu- 
lence capable of mixing individual SST exhaust plumes with . the 
ambient atmosphere appears to depend upon occasional breaking 
waves arising from either external perturbations, e.g., mountain 
waves, or unstable internal gravity waves. As a result, it can be 
assumed that small-scale turbulence in the lower stratosphere is 
both sporadic, rather than continuous, and of relatively short 
duration. Aircraft observations (ref. 3) suggest detectable levels 
of small-scale turbulence occur only one to five per cent of the 
time, and individual turbulence "bursts" persist for only a few 
minutes to an hour or so. 
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If this general picture of turbulence, on the scale necessary 
to diffuse the initial SST exhaust plume, is approximately correct, 
we should expect that a very large fraction of these plumes will be 
emitted into an atmosphere which is in essentially laminar flow 
conditions. However, there is one important source of small-scale 
turbulence which must be considered, namely, the turbulence induced 
by the passage of the aircraft itself, particularly the decay of 
the wing-tip vortices. The sequence of events which this turbu- 
lence produces is shown in Fig. 1. Unfortunately, the quantitative 
treatment of the vortex roll-up, the onset of instability, and 
the structure and decay of small-scale turbulence during the five 
to ten minutes it takes for these events to occur is not yet avail- 
able. In this part of the problem we must still settle for largely 
qualitative and empirical estimates. The importance of this simple 
fact, and the need to acquire a considerably improved resolution 
of this phase of the problem, resides in the chemistry which goes 
on during this brief but critical time, and in the specification 
of the concentration distributions of critical exhaust components 
as they enter the post-vortex phase of their existence (Fig. 1). 

From the viewpoint of the subsequent effects of N0 X and H p 0 
on the environmental O3, we are particularly interested in the 
entrainment of 0^ and the conversion of NO to NO2 during the 

vortex decay phase. During the vortex phase, the SST exhaust 
plume is characterized by concentrations of NO and H^O which are 

orders of magnitude above ambient values. However, during the 
vortex roll-up phase these plumes are also essentially 0-^- and N0 p - 
free. O3, of course, is relatively abundant in the environment 

surrounding the plume, but, until vortex breakup occurs, there is 
no guaranteed source of turbulence to mix the NO and 0 Q . If this 
mixing is initiated and sustained for a few minutes, a^highly 
significant chain of chemical reactions is initiated: 

1. 0^ reacts with MO to produce NOg and 0 2 . 

2. In the sunlit sky, photodissociation of N0 2 and 0 P 
produces 0{3P) and 0('D), as well as NO and Og. J 

3. 0(3 P ) reacts with 0 2 and M to form 0 p . 

4. 0( ’D) reacts with HgO to form 2 OH radicals. 

5. OH reacts with N0 2 to form HNO^. 

Other reactions occur as well, but this sequence, occurring 
as it does . at greatly elevated concentrations of NO and H P 0, poses 
the possibility that a significant fraction of the ^NO is converted 
to HNOq during the first few hours of the plume's existence. Given 
the inert character of HNO^, so far as 0^ is concerned, and its 

slow photodissociation back to NOg and OH, this conversion, if 
realized, strongly buffers the N0 y catalytic cycle for the destruc- 
tion of O3. 
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Figure 1. A qualitative description of the early history of 

individual SST engine exhaust plumes, the turbulence 
fields, and the critical chemical species involved 
in determining the fate of 0 Q and N0 V during this 
time. 5 x 
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When and if this is the case is a major problem considered 
here. For now, however, it is important to realize that this 
whole sequence depends critically upon the conversion of NO to 
NOg in the plume while the concentrations of N0 X and H 2 0 are 

still large . Th i s sequence must be Initiated by the entrainmen t 
of environmental 0^, a process which necessarily leads to the 

dilution of the plume -bound N0 X and H 2 0. Clearly, we have a 
potentially fine balance between diffusion and kinetic processes 
in a system which is characterized by inhomogeneous mixing. 

Against this background, the Aeronautical Research Associates 
of Princeton (A.R.A.P. j has undertaken an evolutionary program of 
model development and analyses in support of the CIAP objectives. 
The initial problems addressed were directed primarily to the rate 
at which individual SST plumes would spread until they achieved 
lateral dimensions of the order of the large-scale (global) models' 
grid spacing. (During this time global models cannot resolve 
individual plume features.) Coupled with this diffusion study, we 
have examined the problem of relatively simple chemistry between 
exhaust and environmental constituents, and to what degree these 
reactions might be limited by either inadequate turbulent mixing 
or by inhomogeneous mixing. The results of these early studies 
have been reported (refs. 2 and 4) and are not repeated here. 

In the latter stages of this study we have focussed attention 
on the problems associated with combined chemistry and diffusion in 
individual SST plumes during the time they exist in high concentra- 
tions and in relatively intact plume segments. Attention has been 
directed first to the problems of interactions between the oxides 
of nitrogen, NO and N0 2 , and 0^ and the impact of turbulent mixing 
on these processes. The results of this study have been reported 
by Hilst (ref. 5). 

As a result of the N0 x -0^ study, it was concluded that the 
possibility of a significant conversion of NO and N0 2 to HNO^ 
during the early plume history deserved a detailed analysis. This 
problem constitutes the final phase of this study for CIAP and, 
since it subsumes all of the earlier analyses, it is the primary 
topic of this report. In addressing this latter problem, it is, 
of course, necessary to include the water vapor component of the 
SST exhaust. This makes it possible to address the important 
problem of the comparative roles of H2O and N0 X in determining the 
0^ balance in the stratosphere. 
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THE ANALYTICAL TOOLS 


Based on the earlier work of Donaldson (refs .6, 7) an ^ 
Donaldson and Hilst (ref. 8), it was evident that the techniques 
of second-order closure of the equations for turbulent dynamics 
and chemical kinetics, using the principles of invariance postu- 
lated by Donaldson, were particularly applicable to the problems 
of early SST plume behavior. In fact, the problem of inhomogeneous 
mixing and its effect on fast chemical reaction rates required an 
ability to handle the prediction of the dynamic behavior of "second - 
order correlations of concentration fluctuations in turbulent mix- 
ing. These correlations could be expected to be strongly time 
dependent, and not subject to simple parameterization by mean value 
characteristics of the flow, plume shape, and chemical reaction 
rates. An analysis of the possible limits of the effects of 
incomplete mixing by Hilst (ref. 2) showed this could be a signi- 
ficant feature of early SST plume behavior, but deferred the 
identification of when or under what circumstances this is the 
case to full model calculations. 

Beyond this requirement for dynamic models which explicitly 
include the second -order moments, it was also recognized that the 
invariant models provide a ready check on the circumstances under 
which first-order closure models provide an adequate representa- 
tion of the combined chemistry and diffusion processes. These 
model predictions, and the parameters associated with them, such 
as effective diffusivity, are recovered explicitly as limiting 
cases of the second -order closure model predictions (cf. Donaldson 
(ref. 1)). The exercising of the invariant model for the SST 
plume conditions also served the useful purpose of "tuning up" the 
first-order closure models used by others. These results have 
been summarized by Hilst (ref. 9). 

The detailed development of invariant models for the predic- 
tion of turbulence and turbulent diffusion in shear flows has been 
reported by Donaldson (ref. l); the development of an invariant 
model for the chemical kinetics of binary, isothermal reactions 
and the coupling of these kinetic models with the diffusion models 
have been reported by Hilst (refs. 2, 5). Briefly, these modeling 
techniques, which date back to Keller and Friedman (ref. 10) but 
were implicit in the work of Reynolds (ref. 11 ), take advantage 
of the fact that the dynamical equation for the behavior of two 
jointly fluctuating variables, say x' and y f , can be written 


d x 'y 1 
ciT 


x' 


dy 1 
<5t 


+ y 


dx 1 

3F" 


(l) 
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For example, and as an illustration of the second -order closure 
technique, consider the chemical kinetics of an isothermal binary 
reaction specified by 



k 

a + (3 -*• y + 6 

(2) 

Then 


= - k[a][P] 

(3) 

and 

||£1 =■ - KCa] [p] 

w 


where [i] denotes the local concentration of the ith chemical 
species, and k is the appropriate reaction rate constant. For 
convenience of notation let x = [a] and y = [P], and we assume 
ensemble averages of x and y can be formed such that 


x = x + x 


and 

where jP" = yl" = 0, 
But 


y = y + y* 


dx 1 dx 8x 


3F _ 3F 


From Eqs. (3) -(6) 


§21 _ 

" 3F cfiT 


8x <3x 
3t + St~ 


- - k(x y + X y' + y x' + x'y') 


and 


5y 8y 

3t + 3^ 


= - k(x y + x y» + y x' + x'y 1 ) 


Taking the ensemble averages of Eqs. (9) and (10) 


57 

3t 


= - k(x y + x 'y ' ) 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


( 11 ) 
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and 


( 12 ) 


= - k(x y + x'y* ) 
Subtracting (ll) from ( 9 ) and ( 12 ) from ( 10 ) gives 




dx ' 

ST" = " 

k(x y' + 

y 

x ’ + x 'y ' 

- x'y 

f ) 

( 13 ) 

and 


dy ' 

TF* = " 









k(x y' + 

y 

x ' + x'y' 

-x'y 


( 14 ) 

Multiplying 

(13) and 

( 14 ) by 

y ' and x 

1 

7 

respectively, 




y ' 

dx ' 

= ' 

k(x y ' 2 - 

(- y x 1 y ' + x 

lyl 2 _ 

y 'x'y 1 

r ) 









( 15 ) 

and 










x ' 

dy 1 

‘ 

k(x x'y 1 

+ 

- .2 
y x ' + x 1 

| 2 y» _ 

x'x 'y 1 

r 


(16) 


Adding (15) and (16) and again taking an ensemble average, 

= - k(x y' 2 + y x »y 1 + x 1 y 1 


- k(x x'y 1 + y x 1 2 + x' 2 y') (17) 

Equation ( 17 ) introduces four new variables, x' 2 , y' 2 , x ' 2 y ' , 
2 

and x'y 1 , for which prediction equations are required. The first 
two are readily obtained by multiplying (13) and (l 4 ) by 2x 1 and 
2y ' respectively, and averaging 


2k(x x'y' + y x ' 2 + x 1 2 y ' ) (18) 


2k(x y ' 2 + y x'y' + x'y' 2 ) (19) 


and 


dx ' 1 
3F~ 


dy '' 
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At this stage, the chemical kinetics model is 


and 


dx 

3t 53 



dx l y 

cTE 


dx 1 ' 

= 



- k(x y + x'y 1 ) 

- k(x y + x'y 1 ) 


Q p 

- k(x y 1 + y x 'y 1 + x 'y 1 

k(x x 'y 1 + y x 1 + x 1 y') 


- 2k(x x'y' + y x' 2 + x'^y 1 ) 


- 2k(x y' 2 + y x 'y 1 + x'y' 2 ) 


(20) 

( 21 ) 


( 22 ) 

(23) 

(24) 


a set of five simultaneous partial differential equations involv- 
ing seven unknowns. If we choose to "close u the system at this 

p O 

level, we are obliged to find expr essions for x* y' and x’y 1 in 

terms of x, y, x'y', x' 2 , and y' 2 . Since the first and second 
moments of x and y cannot generally contain all of the information 

2 2 

necessary to specify x' y 1 and x’y' , these expressions are 
necessarily approximate, a situation which is generally true at 
any level of closure. 

We seek these approximate relationships from the general 
definition of a third-order moment, i.e., 

~2— 1 N - ? _ 

x ' y' = n 2 n i^ x + x l) (y + y{) (25) 

Expanding the right-hand side of (25 ) , rearranging terms, and 
dividing by x 2 y gives 
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1+2 


x y \ 


— r-r 72 

x *y 1 , x ' 

_____ + _ 2 

x y x 


-2 

x 


1 ^ 

y N 1=1 


n.x.y . 

1 i J i 


(26) 


p 

(The definition of x'y 1 is obtained by simply interchanging x and 
y in (26).) We note immediately from Eq. (26) that the explicit 
information regarding the third -order moments contained in the 
first- and second-order moments can be identified by the bracketed 
term on the right-hand side, and the problem is reduced to finding 
a useful approximation for the summation term in Eq. (26). This 
term is, of course, the third -order moment of the joint frequency 
distribution n 1 (x i ,y i ) about the origin of the (x,y) plane. 

Prom the general properties of n i (xj_, y^), we can immediately 
deduce the following conditions which Eq. (26) must satisfy: 

1. When either x' or y 1 =0, 


_i_ I 

x 2 y N i=l 


n i x i y i 


= 1 


2. When n^ is nonzero only for (0, y^ and (x^, 0), 


3f 2 


y N i=l 


n i7i 


= 0 and x'y' 


x y 


Then 



1 


N 

2 


y N i=l 


V? y i 


? 

x» y 1 

-2 - 

x y 



1 


(This is the condition for termination of the reaction due to local 
depletion of either of the reacting chemical species.) 

3. When x’y 1 = 0 (random distributions of x^ and ) 
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N 


x 2 y N i=l 


2 n.x.y. 
111 


x ,2 y ' 

-2 - ' -2 
x y x 


x 1 n 
+ ^ — + 1 


4. When x'y' = 0 


N 


1 n 0 x 'y ' x 1 

_ — 2 n.x.y =1+2 _ + ^p- 

x y N 1=1 111 x y x 


The se fou r conditions should be satisfied by any approximate model 
of x |2 y ' . 

2 2 

A useful model of x 1 y' and x'y' , developed out of these 
considerations at A.R.A.P. (ref. 2), is 


-2 - r 


x« 2 yi - 21 ... y . . 

y _ M + 1 


1 + + 2 


tv i 


-2 

x 


x y 


x 'y 
x y 


- M 


and 


- — y_.. „ 

y - M + 1 ^ 


! 1 + ^T5— + 2 


X y 


x'y 1 
x y 


- M 


(27) 


(28) 


where M is either 0 or 1 according to whether 

zzVi 

x 2 y 2 )> 


(29) 


The importance of retaining the higher-order moments in the chemical 
kinetics model is demonstrated for the case when x^ and y^ are 

jointly log-normally distributed with an initial value of 

x' 2 y ' 2 /x 2 y 2 = 8100. The predicted reaction rates , dx/dt, for 
three cases^ 1) mean -va lue chemistry, i.e., x 1 = y 1 = 0, 2) when 

only x'y 1 , x' 2 , and y |2 are retained, i.e., x ,2 y' = x'y' 2 = 0, and 
3) the full A.R.A.P. chemical model, are compared with the exact 
solution of the chemical kinetics equations in Pig. 2. The neglect 
of all fluctuations grossly underestimates the reaction rate; 
neglect of the third-order moments results in a large overestimate 
of this rate; and inclusion of these terms as modeled here captures 
almost exactly the time history of the react ion rate fo r th is 

situation. It should also be noted that as x ,2 /x 2 and y' 2 /y 2 
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Kt 


Figure 2. Comparison of predicted vs actual chemical 
reaction rates under the following condi- 
tions : 

1) Neglect of concentration fluctuations 
(uniformly mixed). 

2) Neglect of third-order correlations. 

3) Inclusion of approximate estimates of 
third-order correlations. 
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become less than 1, the fluctuation terms in the chemical kinetics 
equations become much less than the mean value terms and we recover 
the classical mean-value chemical kinetics model as a limiting 
condition on this more general model (ref. 2). 


The importance, or lack of importance of these fluctuation 
terms will be apparent after the chemical submodel is coupled with 
the second-order turbulent diffusion model. We may . note here, 
however, that the primary effect of inho mogeneous mixing is propa- 
gated through the ’’mixedness" term, x'y 1 . This correlation arises 
primarily through the "folding" of one reactant into the other, 
without the local, intimate mixing necessary for them to react 
chemically. When intimate mixing has been achieved, this effect 
disappears; it has been o ne of our purposes to assess both the 
magnitude and duration of x'y 1 in a realistic turbulence and plume 
configuration situation. 


THE MODEL DEVELOPMENT 


The Essential Chemistry of the Highly Perturbed 
Hydrogen-Oxygen-Nitrogen Stratosphere 


Recalling that our primary purposes in the present study are 
to focus attention on the roles of N0 X and H 2 0 in the determina- 
tion of the balance of 0^, and to examine the fate of N0 X during 
the first few hours after N0 X and HgO are emitted by a single 

aircraft, we have chosen the following chemical reactions as the 
essential chemistry of this highly perturbed region of the lower 
stratosphere : 

k 2 

(Rl) NO + 0^ -*■ N0 2 + 0 2 


( R2 ) N0 2 + 0 -► NO + 0 2 

J 2 

( R3 ) NO 2 + h v — ► NO + 0 


k„ - 4 x 10 


cnr* 


molec - sec 


-12 3 * 

k = 9.1 x 10 crn J 

-> molec - sec 


j = 0 to 6.4 x 10 0 sec 


(R4) 


0 + o 2 



+ M 


k 1 [0 2 ][M] = 2.3 x 10 2 sec' 1 


l4 



** 


1 -I 

o ± 

( R5 ) 0^ + hv — 0 + 

Jl* 

(R5 ') 0^ + hv 0( 'D) + 0 2 


(R6) 0( >D) + H 2 0 - 2 OH 

**■ 6 

(R7) 0( 'D) + M -*■ 0 + M 

Kli 

(R8) OH 4 NO 2 - HN0 3 


(R9) OH + 0 3 - HO 2 + 0 2 

J 3 

(RIO) HNO^ 4 hv - OH + NO 2 

k 8 

(Rll) OH 4 HNO^ - HpO + NO ^ 
(R12) HO 2 4 0, - OH 4 2 O p 


J 1 = 0 to 4 x 10’3 sec” 1 


k„ = 


kr = 


k t , = 


,-10 


,-11 


■12 


*7 


= 2.7 x 10 


-14 


,-13 


k g = 3.2 x 10 


-15 


-4 -1 

sec 

** 

3 

crrr 

* 

molec - 

sec 

3 

cm J 

* 

molec - 

sec 

3 

cm J 

* 

molec - 

sec 

3 

cm J 

* 

molec - 

sec 

10~^ sec 

-i * 

3 

cnr 

* 

molec - 

sec 

3 

cnr 

* 

molec - 

sec 


This reaction set for twelve chemical species represents, of 
course, a rather severe truncation of the complete 0^ chemistry in 

the stratosphere (ref. 12), and the validity of this approximation 
has not yet been adequately tested against the more complete 
reaction/chemical-species simulation models. Until this is done, 
the results derived from the present chemical submodel must be 
properly qualified. 

*Garvin (ref. 13) 

'X' ■X - 

Crutzen (ref. 14) 
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However, these choices of essential chemistry agree well with 
choices made by other CIAP modelers, extending Lockheed's chemistry 
model (ref. 15) , and catching the essence of the LLL model (ref. 16) 
for the plume phase chemistry. 


The Chemical Submodel 

Assuming then that the above reaction set does capture the 
essential chemistry in a strongly perturbed hydrogen-oxygen- 
nitrogen stratosphere, we may write the rate equations for each 
chemical species as follows: 

d[CU 

3^- = k 1 [0][0 2 ][M] - k 2 [N0][0 3 ] - (J x + J 1 *)[0 3 ] 


- k ? [0H][0 3 ] - k 9 [H0 2 ][0 3 ] (30) 


1^-= J 2 [N0 2 ] + J 1 [0 3 ] + k 6 [0(«D)}[M] - k 3 [N0 2 ][0] 

- k 1 [0][0 2 ][M] 

J 2 CN0 2 ] + k 3 [N0 2 ][0] - k 2 [0 3 ][N0] 

^[NO ] 

3 t ' = J 3 [ HN0 3 ] + k 2 [H0][0 3 ] - j 2 [N 0 2 ] - k 3 [N0 2 ][0] 

- k 4 [N0 2 ][0Hj 

d[HN0_] 

^ = k 4 [N0 2 ][0H] - j^HNO ] - kg[HN0 3 ][0H] 


(31) 

(32) 


(33) 

(34) 
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= 2k^[0( «D)][H 2 0] + J 3 [HN0 3 ] - k 4 [N0 2 ][0H] - k ? [0H][0 3 ] 


+ k 9 [H0 2 ][0 3 ] - kg[HN0 3 ][0H] 


( 35 ) 


b[ 

37 


0(«D)1 


J 1 *[0 3 ] - k 5 [H 2 0][0('D)] - k 6 [0( -D) ] [M] 


(36) 


5[h 2 0] 

"St 

S[NO ] 

S[o 2 ] 

"St 


S[ho 2 ] 
3t 


kg[HN0 3 ][0H] - k 5 [0{'D)][H 2 0) 
kg[0H][HN0 3 ] 

k 2 [ NO ] [ Og ] + k 3 [N0 2 ][0] + k 7 [0H][0 3 ] + 2k 9 [H0 2 ] [ 0 3 ] 

} 

(J x + J!*)[0 3 ] - fc L [0][0 2 ][M] 

k ? [0H][0 3 ] - k g [H0 2 ][0 3 ] 



(37) 

(38) 


(39) 


(40) 


(4i) 


This set of twelve simultaneous partial differential equations 
represents the nonlinear, closed set which can be derived from 
reactions (El) through (R12). This full set can be solved by 
numerical techniques, given appropriate input values for the twelve 
chemical species concentrations and the values of the reaction rate 
constants. Some of the latter are time and position dependent, 
particularly the photodissociation "constants," but we choose to 
suppress this variability for our initial calculations. This choice 
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limits our time integration period to, at most, a few hours, and 
either sunlit or dark skies. 

Before proceeding to these calculations, however, we may take 
advantage of several approximations which are appropriate to the 
SST near wake and which greatly simplify both the chemical calcula- 
tions and the later coupling with source, diffusion, and sink terms. 
Since some of these approximations constrain the regions of validity 
of the chemical submodel, they should be carefully observed. 

First, we note that the near wake plume is expected to be rich 
in water vapor and molecular oxygen. We therefore set Eqs . (37) 
and ( 39 ) equal to zero, and treat [H^O] and [Og] as constants in 
the chemical model. 

Turning attention now to reactions (r 8), (R9), (Rll) and (R12), 
we note that, so long as [N0 2 ] > 0.01 [0^], the completion for OH 

favors the production of HNO^ over the production of H0 2 . Later, 
when a substantial amount of N0 2 has been converted to HNO 3 , the 
preferred reaction for OH is with HNO^, to form H 2 0 and N0^ 

(reaction (Rll)). This reaction effectively removes some of the NO 
and N0 2 from the system. On this basis, we suppress reaction (R12) 

as a negligible part of this complex system and delete the term 
k g [H0 2 ][0 3 ] in Eqs. ( 30 ) and (35). 

For the next set of approximations, we note that 0 and 0( >D) 
are highly reactive chemical species and can be expected to be in 
approximate equilibrium at all times. Assuming this is the case, 
we may set Eqs. (31) and ( 36 ) equal to zero and solve for the 
equilibrium concentrations of 0 and 0('D). This yields 


[0] 


J 2 [M0 g ] + 3 1 [0-^1 + k 6 [°( '!>)][«] 

kjJOgKM] + kjINOg] . ' 


(42) 


and 


[0( 1 L) ] = 


j-i*[0 Q ] 


k 5 [H 2 0] 


kr[M] 


(43) 


But the maximum value of [H 2 0] expected in the post vortex plume 
is about 10 1 5 to lO 1 ^ molecules/cnP and, taking [M] = 10 mole- 
cule s/cm3, k 5 [H 2 0] « kg[M] . Then 


[0( 1 D) ] - 


Jl* [0 3 ] 


^[M] 


^This approximation is valid for [H 2 0] < 10^° molec/cnv 5 


( 44 ) 


18 



Substituting Eq. (44) in Eq. (42), and noting that - 10 ^ , 


[0] 


jg[N0 2 ] + J 1 [0 3 ] 

fcJOgHM] + k 3 [N0 2 ] 


( 45 ) 


Equation (45) can be simplified further by noting that, for. 
[NOg] < 10 13 (the range of interest), ‘k^NOg] < k 1 [0 2 ][M], and 


[0] 


jg[NQg] + J I [Qg3 

k 1 [0 2 ][M] 


(46) 


In arriving at these approximations the primary constraint is 

.01[0 3 ] < [NOg] < 10 13 molec/cm 3 

a range of values of [NOg] which must be observed in the calcula- 
tions . 

Using all of the' above approximations, Eqs. (30) through (4l) 
may be written: 

j 2 [N0 2 ] - k 2 [N0][0 3 ] - k- [OH] [0^1 ( 4 ?) 


j 2 [N0 2 ] + k 3 [N0 2 ] 


j 2 [nq 2 ] + j x [o 3 ; 

k x [0 2 ][M] 


- k 2 [0 3 ][N0] (48) 


*[0 3 ] 

St 


s[no; 

st 


s[no 2 ] 

St 


j [HNO^ + k 2 [N0][0 3 ] - jg[N0 2 ] 


k 3 [N0 2 ] 


j 2 [N0 2 ] + J 1 [0 3 ] 
k x [0 2 ][M] 


- k 4 [N0 2 ][0H] 


(49) 
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(50) 


d[HN0,] 
3t * L_ 


k 4 [N0 2 ][0H] - j 3 [HN0 3 J - k g [HN0 3 ][0H] 


a 


OH] 


2k j * 

[0 o ][H 0] + jjHNO ] - kjNO ][0H] 
k 6 [M] 3 2 3 3 42 


- k 7 [0 3 ][OH] k g [ HNOg ] [OH] 


(51) 


Equations (47) through ( 51 ) have been solved for the variety 
of initial conditions which might be expected in the post-vortex 
SST plume for a period of one to three hours. These calculations 
(which are appropriate to nonturbulent conditions in the strato- 
sphere) are discussed in the next section. We now complete the 
model construction by coupling the chemistry submodel with the 
diffusion model . 


The Coupled One -Dimensional Diffusion/Chemistry Model 


We adopt a frame of reference which moves with a volume of 
the plume segment and note that mass continuity for the jth chemi- 
cal species, requires that 



5 g rn 

3x? 


U ] [ j ] 


5 

3t 


in) 


chem 


+ S. 


(52) 


where the x^ denote orthogonal directions, u| are the components of 
turbulent motion along these directions, and Sj is any local source 
or sink for j, other than chemical alteration. We take Sj = 0 for 
our present problem, and confine attention to turbulent diffusion 
in the lateral (y) direction only. Then, 


d[ j] 

dt 



v ' [ j ] ' + 


5 


[j]) 


chem 


(53) 


and our problem reduces to finding appropriate models for each term 
on the right-hand side of Eq. (53) and for each of the five chemi- 
cal species retained as variables in the chemical submodel. These 
modeling techniques are provided by previous development of second- 
order closure techniques for turbulent diffusion processes (ref. l) 
and inhomogeneous chemistry (ref. 2). 
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However, > in order to keep the number of simultaneous equations 
manageable, we invoke two more approximations in the chemistry 
model. First we assume that OH is always in equilibrium, and 
second, we neglect the second-order correlations of concentration 
fluctuations in the chemistry, i.e., we use mean value chemistry 
only. The importance of this latter assumption can be rigorously 
tested only when it is deemed desirable to program the more complete 
model for a large, high speed computer. However, calculations on 
simpler reaction systems (HoO is neglected) in which the second- 
order correlation, or "mixedness" term is retained have shown this 
term alters the reaction rate predicted by mean value chemistry by 
no more than about thirty per cent for typical SST plumes and, in 
turbulent flow, persists for only a matter of a few tens of minutes 
(see Appendix A). Given the degree of uncertainty in other parts 
of the system, e.g.g., reaction rate constants, initial species 
concentrations, turbulence intensity, etc., we are justified in 
neglecting this effect in these preliminary calculations for the 
full H-N-0 system. 

The justification for assuming that OH is always in chemical 
equilibrium must be based on its high degree of reactivity and 
therefore its small characteristic chemical time in comparison 
with 0-,, NO, N0 2 , and HNOg. Since this assumption is not made in 
the full chemical submodel, we are in a position to assess its 
validity by comparing calculations made with that model with those 
made with the coupled model when the turbulent diffusion is taken 
as zero. Adopting this assumption, we have from Eq. (22), 


2 k j * 

k-flY[] + j^HNOg] 

[OH] = -2 - — (54) 

k 4 [N0 2 ] + k ? [0 3 ] + kg[HN0 3 ] 

A comparison of k^, ky, and kg shows that, so long as [N0 2 ] > .01[0g] 

or [ HNOg] > 0.1[0 3 ], we may safely neglect the term ky[0g] in 

Eq. (54), although we retain it in Eq. (47). Finally, we note that 
in the early plume phase, where [OH] can be expected to be relative- 
ly large, the primary mechanism for the destruction of HNO^ is by 

reaction with OH, not by photodissociation. We therefore neglect 
j 3 [ HNO g ] in favor of kg[HN0g][0H] wherever they appear. 

Under these conditions, the coupled chemistry /diffusion model 
becomes : 
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d[o ] & 

= " By v ^ 0 3 ]' + J 2 [N0 2 ] - k 2 [N0][0 3 ] 



2k 5 J 1 *[0 3 ] [H g 0] 
k 6 [M](k 4 [N0 2 ] + kg[HN0 3 ]) 


(55) 


df NO] a 

dt By 


v i [NO] ' 


J 2 [no 2 ] 


^[mo 2 ] 




J 2 [no 2 ; 


J 1 [o 3 ] 


] 


k 2 [0 3 ][NO] (56) 


d[HNO„] 

dt 


By V ' ^ HN0 3 3 


k 4[NQ g ] - kg [ HNO^ ] 
k 4 [N0 2 ] + k Q [ HNOg ] 


2k 5 J 1 *[Q^] [H^Q3 
k 6 [M] 


(57) 


d[NO x ] 

35 


a 

"5y 


v[NO x ] 


- k 


8 


[hno 3 ] 


2k,. J*[Q^3 [H 2 0] 
k: 4 [N0 2 ] + kg[HN0 3 ] 


(58) 


d[H 2 0] v 

d5 = " By v, t H 2 0 ^ ' (59) 


[N0 2 ] = [NO x ] - [NO] - [HN0 3 ] 


(60) 
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and it is understood that all unprimed chemical species concentra- 
tions are ensemble averages. 

Having chosen to suppress the higher-order moments in the 
chemistry submodel, we may write down the diffusion equations 
directly from Donaldson (ref. 1) in the form 


dv i[ JTL 
3E 



+ 7^- Tt[ j] I 



3v ' 

w~ 


JL 

At 


V '[ J] 


(61) 


ST' 


dt 




= v . [ j ] . - v.T. 


IF 




(V* 


stttT 


w 


(62) 


dw i T i 




+ 2 


d / £v <Tt \ 

37 ( ay ) 



-3- v iT t 
*1 


and 


dT» 

3t~~ = ' 


2v iT i 


dT & 
3y + 3y 




) 


(63) 


(64) 


where molecular diffusion has been neglected. Since the chemical 
species do not appear in Eqs, (63) and (64), these are derived 
from the turbulence model and used as standard inputs to the five 
pairs of diffusion equations denoted by Eqs. (6l) and (62). The 
latter are, in turn, coupled with the appropriate chemistry terms 
to provide a set of fifteen differential equations which are to be 
solved for the values of [O3], [NO], [NOg]# [HNO^], and [H2O] as a 

function of position in the plume, y, and time after vortex decay, 
t. The inputs required for the turbulence part of this model a re 

p P 2 

the turbulent kinetic energy, q , the lateral component of q , v< , 
and the scale lengths, A^. The reader interested in the details 

of this part of the model is referred to ref. 1 and 2. The main 
point to note here is that the diffusion scale length is keyed to 
the lateral dimension of the plume, giving a scale-dependent 
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diffusion process analogous to that developed by Walton (ref. 17). 

In the present use of this coupled model we ^re more 
interested in exploring the range of possible turbulence and 
chemistry processes in the near SST wake than in duplicating any 
single event. The lack of validation data is not a severe handi- 
cap in this initial sensitivity analysis; choosing which of the 
many possible conditions actually does occur, and validating the 
model's predictions for those cases, does of course depend 
heavily on data which are not yet available. 


PRELIMINARY MODEL CALCULATIONS 


Chemistry Only 

Since one of the primary problems addressed by CIAP is the 
chemical reactions among N0 X , H 2 0, and O3, particularly as these 
effect the balance of 0^ in the stratosphere, we have exercised 

the chemistry model extensively to determine the rates and time- 
integrated reactions for these components during their first few 
hours in the sunlit stratosphere. Recall one of our primary 
purposes here is to estimate the chemical form of the oxides of 
nitrogen and how they are transformed during these first few hours. 
For these preliminary calculations we have programmed Eqs. (47) - 
(51) for computer solutions and have calculated the time-dependent 
values of [O3], [NO], [NOg], and [HNO^] for all combinations of the 

following initial conditions: 

[0 3 ] o = 5 x 10 12 , 10 12 , and 10 11 molecules/cm 3 

[H 2 0 ] q = 3.5 x 10 17 , 1.7 x 10 17 , 3.5 x 10 l6 , 3.5 x 10 l4 , 
and 3.5 x 10 13 molecules/cm 3 

[N0 x ] q = [N0] q + [N 0 2 ] q = 10 13 , 5 x 10 12 , 2 x 10 12 , 10 12 , 

ll" o 

and 10 molecules/crrr 


[no] o /[no 2 ] q 


0, 1, and 3 


[hno 3 ] o = 0 
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The reaction rate constants were assigned values given 
reactions (Rl)-(R12), and the photolysis constants, j-^. 


-4 


were assigned values of 4 x 10 - ^, 6.4 x 10'^, and 4 x 


with 

<3 2 ’ 


and 


10 ' 5 , 


respectively. Approximately 200 calculations were performed for 
a time integration period of 2.5 hours. Typical time histories 
of chemical species concentrations are shown in Fig. 3J in 
addition to these primary outputs, the program also provided the 
time rate of change of species concentration and the character- 
istic time for these changes, also as a function of time. 


As an initial framework for examining this multi-variate 
system, and for an early appraisal of HN0 o formation in this 

photochemical simulation, we have determined the combinations of 
[O^Jq, [N0 X ] , an ^ [NO] 0 /[NOo] q for which more than ten per cent 
of^the initial N0 X is converted to HNO^in 2.5 hours. These 

delineations are displayed in Fig. 4. The following qualitative 
picture of the probability of significant HNOo production is 
immediately evident: 


1. The probability of HNO-^ formation increases with increas- 
ing [HoO] and increasing [0p] n “’and decreases with increasing 
[N0 X ]. 


2. HNO^ formation is less sensitive to [N0] o /[N0 2 ] o than to 
[0 3 ] o , but the sensitivity to the former increases as [N0 X ] and 

[H 2 0] increase and [0^] decreases. 

3, For [N Oyl > ~ 5 x 10 lg molecules/cm? the probability of 
HNO g formation becomes vanishingly small . 

The combinations of [0,] Q , [N0 X ], and [N0] Q /[N02] o used in 

constructing Fig. 4 bracket essentially all conceivable values 
which could be encountered in the SST exhaust plume. However, as 
a preliminary test of the sensitivity of HNO q formation to the 

photolysis rates, approximately 100 additional calculations were 
made for the same initial conditions, but was increased to 

1.6 x 10*’^ sec -1 and j 2 was decreased to 6.4 x 10 - ^ sec - ''". The 
combinations of [0 3 ] o , [N0 X ], [H 2 0], and [N0] o /[N0 2 ] o which led 
to the conversion of at least ten per cent of the N0 X to HNO^ in 
2.5 hours are shown in Fig. 5. The same general results described 
above are still evident, but the probability of significant HNO^ 

formation is now considerably enhanced for lower values of [O^Iq 
and moderate to large values of HpO. This result, of course, 
reflects the greater abundance of 0( >D) and consequent formation 
of OH, and then HNO^. 
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Concentrations of designated chemical species 
(molecules/cm 3 ) 





molec/cm 



Figure 4. The threshold values of initial [O^], and [N0]/[N0 2 ] 

which produce more than 10 per cent conversion of N0 X to HWO^ 
in 2.5 hr s'. (j-.* = 4xlO -5 sec -1 and J 2 = 6.4xl0 _i; sec X ) . 


[N0 X ] molec/cm 



J.*lH«0] 


molec 

cm 3 -sec 


Figure 5. The threshold values of initial [0^3, J * [ H^O ] , 

and [N0]/[N0 2 ] which produce more than 10 per cent 
conversion of N0 X to HN0 3 in 2.5 hrs. (j * = 

4xl0 _i+ sec -1 and J 2 = 6.4xl0 - ^ sec -1 ). 
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Even with the elevated value of j^*, the controlling upper 
limit of [N0 X ] is still very evident. An inspection of the calcu- 
lations quickly reveals that the reason for this limit is to be 
found in the abrupt depletion of Oo in the photochemical system 

IP ^ q 

when [N0 X ] exceeds about 10 molecules/crnd . in order to show 

this, we have analyzed the first set of calculations (daily 
average photolysis constants) for the combinations of [N0 X ], [H 2 0], 
and [N 0 ] o /[N 02 ] o for which [O 3 ] is depleted to less than one- 

half its initial value in 2.5 hours. These delineations are shown 
in Fig. 6 , and they show the abrupt demise of O 3 for values of 
[N0 X ] > ' 5 x 10 12 molecules/cm^ and J]_*[H 2 0] <~10- 1 -^ molecules/cm^ 


sec. An interesting interplay of H 2 0 and N 0 X sets in when 
j-L*[H 2 0] > lO 1 ^ molecules/cm3 sec. First, the threshold of [N0 X ] 
for O 3 destruction greater than 50 per cent is raised , and then, 
for J]_* [H 2 0] > ~ 5 x 10 12 molecules/cm^ sec, this measure of the 
destruction of O 3 becomes essentially independent of [N0 X ] . 

Evidently, the first role of H 2 O in the O 3 balance is to inhibit 
the destruction of O 3 by N 0 X , by the simple expedient of convert- 

ing N0 X to HNO 3 , and then at-large values of [H 2 0], H 2 0 alone is 
capable of annihilating O 3 through the OH chain. 

The primary preliminary result to be noted at this stage, 
however, is that significant formation of HNO^ can occur with 

values of [O 3 ], [H^O], and'[N0 x ] which could occur and persist 

for 2.5 hours in the Initial SST exhaust plume. However, given 
the sensitivity of HNO 3 formation to the existence of 0 ^ in the 

plume and the upper bounds on [N0 X ] and [H 2 0] we must inquire 

more closely regarding both the initial concentrations and the 
effect of turbulent mixing before any definitive application of 
these results can be made to the problem of the SST plume. 


Before proceeding to this larger problem, a more quantitative 
analysis of the conversion of N0 X to HNO 3 and the destruction of 
0 _, due to chemical reactions only, is possible and instructive. 
Turning to the destruction rate for O 3 first, we consider the 
characteristic time for 0 o destruction, Tq , as defined by 


T 


0, 


1 d ' [ 0 3 ] \ 

[ 0 3 ] dt / 


(^ 5 ) 


29 



[NO x ]mo!ec/cm 



Figure 6 . The threshold initial values of [N 0 X ], J 1 *[H 2 0], and [N0]/[N0 2 ] 
which produce more than a 50 per cent depletion in [O 3 ] in 
2.5 hrs. (j-,* = -4x10“^ sec -1 and j 0 = S.^xlO 1 ^ sec -1 ), 



and as calculated by the chemical model at t = 2.5 hours. The 
behavior of -Tq as a function of [0 3 ] , [N0 x ], and is 

shown in Pigs. 7 , 8, and 9. The interplay of HgO and N0 X in 
controlling the rate of destruction of 0^ is now very clear. 

Three distinct regimes, as delineated by the solid lines in 

Figs. 7 > 8, and 9 are recognized. First for [N0 X ]/[0^] Q < Ji^HgO]/ 

5 X 10 12 , Tq is independent of [N0 X ] and is inversely proportional 

to j^*[HgO]. In the second regime, Tq increases with increasing 

[N0 X ] until the critical value for N0 X destruction of Oo is 
approached. This critical limit of N0 X depends upon both [0 3 ] o 
and «Ji*[H 2 0] in a somehwat more complicated way than was true for 
the delineation of the boundary between regimes one and two. 


We should also note in passing, although not included in ttjie 
calculations, that at normal background values of [0^] = 5 K 10 , 

[NO ] - 10 10 , and J-^tHgO] - 10 8 , the N0 X limit on 0 3 destruction 
clearly predominates over any I^O effect. Again of peripheral 
interest to the intent of the present study, these simulations 
point out that the characteristic time for 0 3 destruction in near- 
background concentrations of N0 X and 0^ is at least of the order 

of thousands to tens of thousands of hours. If these estimates 
have any validity, they point to the importance of the residence 
time of SST generated N0 X in the stratosphere as the single most 

important parameter in estimating the long-term elfects of these 
emissions on stratospheric 0 3 content. 


Finally, these preliminary and exploratory calculations 
provide a more quantitative analysis of the conversion of N0 X to 

HNOo . Again, we choose to examine this process at 2.5 hours after 

vortex decay and for the conditions of daily average photolysis, 
with other conditions chosen to bracket all reasonable combinations 
of [H 2 0] and [N0 X ] in the initial, post -vortex SST plume. The 

fraction of N0 X converted to HN0 3 during a period of 2,5 hours for 
values of [HgO] and [N0 X ] ranging from approximately ten times 
background to one thousand times background, are shown in Figs. 10,. 
11 and 12. A casual inspection of these results points up 
quickly that combinations of large values of [0 3 ] o and [HgO], 

matched with relatively low values of [N0 X ], produce very large 


conversions of N0 X to HN0 3 in 2.5 hours. The larger question is, 

of course, whether or not these combinations can be produced and 
sustained in the early SST plume. In order to further emphasize 
this very important role of the joint values of [ ^3 1 0 ^ [ N0 X J 
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at t = 2.5 hrs (sec) 


[0 3 ] o - 5XI0 12 mol/cm 

[no] 0 /[no 2 ] 0 =o 

j 2 s 6.4X|0" 3 sec -1 


1.4X10 


1.4X10' 



7X10' 


[N0 X ] molec/cnr 


The characteristic time for the depletion of 
[0 3 ] , at 2.5 hrs, as a function of the initial 

values, of [N0 X ], j-^*[H20], and for [0 3 ] Q = 

5 x 10 12 molecules/em- > . 


Tie characteristic time for the deple 
O^], at 2.5 hrs, as a function of th< 

alues of [ N0 X ] , J 1 *[H 2 0], and for [0 

0 12 molecules/em^ . 






[Oj]o = | ° mol/cm 


[no] 0 /[no 2 ] 0 =o 

- l.4X|0 n ^\ 

j 2 «6.4Xlcr 3 sec~' 

\ 

_ 

A 

L4 X I0 12 

\ 

1.4 X 10*3 

/ 

l 


[NOnJ molec/cm' 


Figure 9. The characteristic time for the depletion of 

[0 3 ], at 2.5 hrs , as a function of the initial 

values of [N0 X ], and for [ 0 ^ 3 0 = 

10 11 molecules/cm^ , 




[HNOj /[no,] at t=2.5 hrs 


[Oder IO M I0 12 5XI0 12 


j I*" [ H 2°] = 1 ^ X| ° 13 

. 2 - ([H 2 0]>~5X|0 16 ) 



[NO x ] mo!ec/cm : 


Figure 10. The fraction of N0 X converted to HNO^ m 

2 5 hrs as a function of the initial [NO^.], 
[o'], and j|[H 2 0] = 1.4 x lO^ molecules/cm 


[HN0 3 ]/[N0x] at t = 2.5 hrs 


i*[He0]«l.4XI0 l: 



The fraction of N0 V converted to HNO- In 
2,5 hrs as a function of the initial [NO ], 

[ 0^ ] f and j*[H 2 0] = 1,4 x 10^ molecules/cm3 
sec . 


Figure 11. 


[HN0 3 ]/[N0 x ] ot t=2.5 hrs 



[H 2 0] In the production of HNO^, from .chemistry alone, these 

calculations have been summarized in still another format in 
Fig. 13. Here the combinations of [0^] o , [N0 X ], and [HgO] which 

can chemically produce the conversion of N0 X to HNO^ greater than 
given fractions of [N0 X ] are displayed. These results are, in turn 
summarized in Fig. 14 to show the combinations of [O^Iq and [N0 X ] 
which, in association with two values of [HgO], can convert more 
than 10 and 50 per cent of the initial N0 X to HNO^ in 2.5 hrs. 

These results, although interesting in their own right, are sober- 
ing when considered in the context of the initial SST plume condi- 
tions. They say, in effect, that from self-regenerative chemistry 
alone (no mixing or dilution) , a large fraction of the N0 X is 

quickly converted to HNO^ only if [N0 X ] and [HgO] two to three 

orders of magnitude above background characterize the post- 
vortex plume and, of course, are maintained for 2 to 3 hours. 

It seems clear enough that reliable estimates of the vortex 
decay period dilution of the initial SST plume, in order to 
establish the initial conditions on [O3], [N0 X ], [N0]/[N0 2 ], and 
[ HgO] , now become critical to the present study. In addition, 

estimates of the dilution or redistribution of these quantities 
during the first few hours after vortex decay are also critical, 
and must be coupled with the chemical system simulation if reliable 
estimates of total N0 X conversion to HNO^ within the plume are to 
be had. 


The Problem of Vortex Decay and the Initial Conditions 
for the Post -Vortex SST Exhaust Plume 

Although second-order closure modeling of the behavior and. 
decay of aircraft-induced vortices and submarine wakes in oceanic 
thermoclines shows promise of quantitative predictions in the 
future, it is generally agreed that detailed predictions of the 
turbulence fields associated with vortex-pair breakup and decay are 
presently beyond the state of the art. In order to lend some 
degree of realism to the sensitivity analyses which follow, we 
have attempted to estimate the entrainment of 0^ and the conver- 
sion of NO to N0 2 during the vortex decay period, using the 

coupled chemistry and diffusion model, and approximating the 
decaying turbulence field by a short-lived, but constant value 
of the turbulent kinetic energy. In addition, we have assumed 
that the dimensions of the turbulence field induced by vortex 
decay are large compared with the width of the exhaust plume 
which is embedded in this turbulence field. These assumptions 
are, at best, only gross approximations of the real situation. 
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[NO x ] molec/cm 3 



Figure l4. Threshold limits of the initial values of [0^] and [N0 X ] on 
the conversion of N0 X to HNO^ in average sunlight and in a 
period of 2.5 hrs . 


However, they do permit an estimate of whether or not turbulent 
mixing, sufficient to entrain environmental 0 ^ in quantities . 

required to convert NO to N0 2 and to maintain a near-environmental 
concentration of 0 ^ in the plume for a few minutes, can be reason- 
ably associated with the vortex breakup and decay. The initial 
turbulent kinetic energy q 2 , due to vortex breakup appears to be 
of the order of one to two m 2 /sec 2 ; only values of q 2 less than 

the order of 0.1 m 2 /sec 2 can be expected to last for 2 to 3 minutes. 
However, while hydrostatic stability may quickly suppress the 
vertical component of this turbulence, the lateral component may 
persist appreciably longer. 

As a first, very approximate analysis of the redistribution 
and chemical alteration of 0^, H 2 0 , and N 0 X , we have assumed that 

the initial vortex breakup entrains environmental air abruptly, so 
that there is an initial content of 0 ^ in the plume, and, following 
this initial "gulp, 1 * the plume is axisymmetric , with the pollutant 
concentrations, [HgO] and [NO], distributed in a Gaussian profile 

with a standard deviation of 10 meters. We further assume that the 
combination of vertical plume collapse and the rapid damping of any 
vertical component of turbulence (both of which associate with 
large hydrostatic stability), produce a "lamina" of plume in the 
vertical dimension and the pollutant concentrations in this dimen- 
sion are more or less uniformly distributed within the vertical 
plume bounds, but retain their essential Gaussian profile in the 
lateral dimensions initially. 

The next step in this analysis is to estimate the total HgO 

and NO emitted by an SST per unit distance of travel. For this 
purpose we have chosen to use a single value of HgO emission, 

based on an emission index of 1.3 and a cruise mode fuel consumption 
rate of 10,000 kg/hr/engine (these estimates seem to be relatively 
firm). However, the uncertainties in the estimates of NO emissions 
(as well as a desire to investigate the effects of raising or' 
lowering NO emissions), have led us to choose four values of the 
NO emission index ranging from 0.010 to .0005. Finally, it was 
assumed the SST was flying in the sunlit sky, and the photolysis 
constants were assigned their average daily values. 

Six levels of turbulence, ranging from zero to v ' 2 = 2 m^/sec^ 
were run in order to investigate this parameter’s influence on 
plume dilution and chemistry. It was assumed that the chosen level 
of turbulence was sustained for 200 seconds. All of these estimates 
and choices led to the nine combinations of input parameters shown 
in Table 1. In order to facilitate comparison with first-order closure 
models, we have included the value of effective eddy diffusivity, 
ky, which associates with the chosen turbulence fields at t 3 200 sec. 

Since, as noted in the explanatory remarks in Table 1, h. is in fact 
a function of time during this stage of plume growth, these values 
of ky should be used very discreetly. 
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TABLE 1 



ky Is estimated from ky = (Oy - Oy ) /2t at t = 200 sec, 

where ap is the calculated variance of the lateral 
distributions of [HgO] and [NO^], both of which are 
conservative chemical species in this model. Oy is, 
of course, more nearly proportional to t^ than to t 


during this period of time, hence ky is in fact a 
function of time. 
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The pertinent results of these simulations of vortex decay- 
mixing and chemistry are listed in Table 2. and displayed in 
Pigs. 15 and 16 . In Fig. 15 we examine the dependence of the 
reactant concentrations on the level of turbulence ([N0 X ] is held 

fixed), and in Fig. 16 we examine the dependence of these plume 

parameters on the amount of N0 V in the plume (v 1 is held fixed). 

We are particularly interested in the amount of NO converted to 
NO 2 J the concentration of 0^ in the plume, and the concentration 

of H^O. Inspection of Fig. 15 shows that, for an initial center- 
line concentration o f N O of 10^3 f more than half of the NO is 

75 o p 

converted to NO 2 if v 1 is greater than 0.02 m /sec , and that the 
0-^ concentrations in the plume are wel l a bove 10- 1 - 2 molecules/cm^ 
for these levels of turbulence. When v' 2 is in the^range of 0.02 
to 0.1 m 2 /sec 2 the H 2 O concentrations are around lO 1 ^ molecules/cm^ 

These results suggest that very modest levels of vortex decay 
turbulence are capable of setting the stage for HNO^ formation in 

the post-vortex phase of the plume life. However, it is also 
clear that very vigorous turbulence will dilute the H 2 O concentra- 
tions to values which are marginal for HNO^ formation. 

Turning attention to the effects of varying the N0 X content 

of the vortex plume, Fig. 16 shows very clearly that a turbulence 

level of 0.1 m 2 /sec 2 is not sufficient to maintain 0, in the plume 

PI ^ 

when the columnar amount of N0 X exceeds about 10 molecules. 

Under these circumstances, the N0 X is increasingly dominated by NO 
(conversion to N 0 2 is minimal), and there is an increasing tendency 
to form an "Og-hole" in the core of the plume. Since these are 

realistic ranges for the initial NO content of the SST plume, we 
are once again reminded that the expected values of N0 X concentra- 
tions are in a very critical range for 0^ destruction, a result 

first suggested by the prior chemistry calculations. However, for 
our present purposes we simply recognize that the N0 X emission 

rate is a critical parameter,, and choose four initial conditions 
for the post-vortex chemistry and dif fus ion calculations. In 

2 p n 

making these choices we have assumed v 1 = 0.1 mysec^ is a reason- 
able effective value during vortex decay, and have extracted the 
pertinent initial conditions for post-vortex calculations from the 

2 2 2 

vortex decay simulations for v' = 0.1 m /sec and [NO x ]j = 

1.3 x 10 20 , 2.5 x 10 20 , 1.2 x 10^, and 2.5 x 10 2 ^ molecules. All 
of the pertinent input parameters for post-vortex calculations are 
shown in Table 3 . These choices are hardly exhaustive of all 
possible conditions; they do tend to bracket reasonable conditions, 
however, and provide at least a preliminary analysis of the 
immediate post-vortex chemistry and diffusion of SST exhaust plumes. 
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[NO^ molecules 

Figure l6a. The calculated initial post-vortex centerline values of [N0 X ], 

[N0 2 ] for a fixed value of turbulence and as a function of the 

NO in a lateral plume column of unit cross-section. 

.X 

Figure l6b . The calculated initial post-vortex columnar amounts of N0 X , NO 

no 2 . 


,22 

[NO], 

total 

, and 


Integrated Concentrations 


V 2 

m 2 /sec 2 

[ NO , 

J 

Jl 

[NOlj 

tnolecules 

[NO 

2^1 

.10 

1.26 

X 

10 20 

4.4 

X 

IO 19 

8.0 

X 

10 19 

0 

2.52 

X 

10 20 

1.5 

X 

10 20 

1.0 

X 

10 20 

.02 

2.52 

X 

10 20 

1.1 

X 

10 20 

1.4 

X 

10 20 

.10 

2.52 

X 

10 20 

8.7 

X 

IO 19 

1.6 

X 

10 20 

.50 

2.52 

X 

10 20 

7.4 

X 

10 19 

1.8 

X 

io 20 

1.00 

2.52 

X 

10 20 

7.1 

X 

10 19 

1.8 

X 

10 20 

2.00 

2.52 

X 

, 10 20 

6.8 

X 

10 19 

1.8 

X 

10 20 

.10 

1.26 

X 

10 21 

7.9 

X 

10 20 

4.7 

X 

io 20 

.10 

2.52 

X 

10 21 

2.0 

X 

10 21 

5.6 

X 

io 20 


2 


'ions and amounts at end of 
reaction. Initially all of 
.tions are listed in Table 1. 


Centerline Concentrations 

t H 2 °>0 

molecules 


[NO ] 

1 x J o 


f°3^o 


1.1 x 10 12 

1.0 x 10 13 

4.3 x io 12 

2.2 x 10 12 

12 

1.0 X lO 1 ^ 

7.6 x 10 11 

5.4 x 10 11 

1.1 x 10 13 

2.2 x 10 13 


7.6 x 10 15 

3.5 x 10 16 

1.5 x 10 16 

7.6 x 10 15 

3.7 x 10 15 

2.6 x 10 15 

1.9 x 10 15 

7.9 x 10 15 

7.7 x 10 15 


12 

3.2 X 10 ^ 

6.7 x 10 11 

2.1 x 10 12 

3.2 x lO 1 ^ 

4.1 x 10 12 

4.3 x 10 12 
4.5 x 10 12 

8.7 x 10 11 

3.2 x 10 11 




TABLE 3 


Initial chemical concentrations and plume dimensions chosen 
for post -vortex plume chemistry and diffusion calculations. 


Case 

Iter/ 12 3 4 

1.88 x io 3 1.88 x 10 3 1.83 x 10 3 1.83 x 10 3 


[0 3 ] q 3.2 x 10 12 
[H 2 0 ] o 7.6 X IO 15 
[N0 x ] q 1.1 x 10 12 
[NO] 3.8 x 10 11 
[N 0 2 ] o 7.0 x 10 11 

[HgOlj 8,8 x IO 23 
[N 0 x ] I 1.3 x 10 20 
[NOJj 4.4 x 10 19 
[N 0 o ] z 8,0 x 10 19 

[hno 3 ] i .9,5 X 10 16 


3.2 X 10 12 8.7 

7.6 x 10 15 7.6 

2.2 x 10 12 1.1 

7.6 x IO 11 7.8 

1.4 x 10 12 3.1 

8.8 x IO 23 8.8 

2.5 x 10 20 1 .'2 

8.7 x 10 19 7.9 

1.6 x 10 20 4.7 

1.9 x IO 13 6.6 


x 10 11 3.3 x IO 11 
x 10 15 7,6 x 10 15 
x IO 13 2.2 x 10 13 
x 10 12 1.9 x IO 13 

Ip 12 

x 10-^ 2.7 X 10 

x IO 23 8.8 x 10 23 
x 10 21 2.5 x 10 2i . 
x 10 20 2.0 X 10 21 
x 10 20 6.0 x IO 19 
x IO 17 1,3 x IO 18 


o 

r i 1 = Plume centerline concentration, molecules/cnv-’ 

L o 

[i] = integrated amount across the plume in a column of unit 

^ cross-section, molecules. 

O 2 

o = Variance of initial [H o 0] and [NO ] distributions, m . 
y ^ A 

The environmental values of [0„] = 5 x 10 12 molecules/cm 3 ; all 

other [i] e =0. ^ 
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COUPLED MODEL CALCULATIONS FOR POST-VORTEX CONDITIONS 


With the. above estimates of the likely range of concentrations 
and concentration distributions for 0^, NO, N0 2 , and HgO at the end 

of the vortex decay period, we may proceed to the primary problem of 
estimating the fate of N0 X and 0^ in the early hours of the post- 

vortex phase. For these calculations we continue the assumption 
that diffusive mixing in the vertical dimension is negligible in 
comparison with lateral diffusion. And, of course, we are inter- 
ested in a range of lateral turbulence intensities, starting with 
no turbulence. In addition, since photolysis is an important 
process in the photochemistry of SST plumes, we have varied the 
values of j|, j 1 , and j 2 for several calculations. 

The combinations of v 1 2 and used in connection with the four 
initial concentration distributions listed in Table 3 are shown in 
Table 4. (j-^ and were varied in proportion to the variation in 

j£.) For these initial analyses we have a total of twenty calcula- 
tions. Each of these calculations was carried out to a real-time 
simulation of 2.5 hours. 

The full computer output for the coupled model comprises an 
overwhelming mass of information on the state of the system and the 
rates at which chemical reactions and diffusive processes are going 
on at each time step for which a print -out is requested. For 
example, a full concentration profile for each chemical species, 
the profile of reaction rates for each chemical reaction term in the 
model, the diffusive flux of each chemical species, and the inte- 
grals and higher-order moments of each of the concentration distri- 
butions are available at each time step. This detail, of course, 
permits microscopic examination of any component of the simulation 
system. Since we are primarily interested in the fate of N0 X and 
0o here, we will forego discussion of these details and concentrate 
on the distributions of 0^, NO, N0 2 , and HNOo at the end of 2.5 hrs 

and as a function of turbulent mixing, photolysis, and the initial 
conditions for reactant amounts and distributions. 

However, as an illustration of the detail with which the 
processes affecting 0^ and N0 X concentrations are portrayed, we 

have extracted from one run (at t = 1000 sec) the profiles of the 
rates at which 0^ and NO are being produced or destroyed by chemis- 
try. These are shown in Fig. 17. Without going into detail, we 
note that in this particular circumstance the conversion of NO to 
N0 2 by reaction with 0^ is completely offset by the photolysis of 

N0 2 and the reaction of N0 2 with 0 at the core of the plume. In 

this location NO is being produced, not destroyed, by chemistry. 
Similar portrayals of the diffusion and other chemical reaction 
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Combinations of initial conditions and turbulence chosen 
for post-vortex plume chemistry and diffusion calculations. 


Turbulence 


Case 


level 1 2 3 4 



m 2 /sec 2 1.6x10 ^ 8x10'^ 1.6xl0~^ 8x10""-* 1.6xl0~^ 8x10”^ 1.6xl0~^ 8xl0 - ^ 


0 xx xx Xx xx 

xx xx XX xx 


x 


.01 

.02 

.05 


x 


X 





Lateral distance from plume 



terms are entirely feasible, if one wishes to pursue such detail. 

As a first general analysis of the fate of 0 ^ and N 0 X , we shall 

confine attention to the concentrations of these materials at the 
centerline of the plume and the integral amount of each chemical 
species in a lateral column of unit cross-section extending through 
the plume. From these measures we can readily deduce the produc- 
tion or depletion of each reactant and the dependence of these 
integral changes on the input conditions. 


The Fields of Turbulence in the Lower Stratosphere 

Before proceeding to these analyses, it is important to 
establish the likely levels of lateral turbulent velocities which 
will operate on the post-vortex plume. Quite literally, the 
natural small-scale lateral turbulence in the lower stratosphere 
is unknown. Inferences can be drawn from aircraft measurements of 
the vertical gust velocities, and, more indirectly, by considera- 
tion of the dispersion of materials in this region of the atmos- 
phere (ref. l8). Neither of these are appropriate to the specifi- 
cation of sustained turbulence on the scale of the post -vortex SST 
plume . 

Given this situation, we have chosen to investigate the 
sensitivity of the fate of 0^ and N0 X to turbulent mixing over a 

range of turbulence intensities which probably bracket any real 
circumstance and which, at any rate, go beyond the value of turbu- 
lence at which chemical alterations are controlled by diffusive- 
mixing. The need for adequate measurements of stratospheric turbu- 
lence, so one can specify where in this range the real processes go 
on, seems obvious, and is highlighted by the analyses which follow. 

As the primary data base for the analyses of the fate of Cb 
and N0 X , we have compiled the values of plume width, as measured 

by the variance of the conservative HgO distribution, <jy, the plume 
centerline concentrations of O3, HgO, N0 X , NO, and NOg, designated 
by [i] 0 > an( ^ the integrated columnar amounts of HgO, N 0 X , NO, NOg, 
and HNO^, designated by [ i ] x , at t = 2.5 hours and for each of the 

twenty runs. These are tabulated in Table 5a - f. In addition, 
the effective eddy difusivities associated with each level of turbu- 
lence used and for t = 2.5 hrs are tabulated in Table 6. Finally 
the calculated deficit of 0^ in the plume, as a fraction of the 

total 0^ which would occupy the plume column if 0^ were at ambient 
concentrations, and the deficit of the 0^ plume centerline concen- 
tration as a fraction of the ambient 0^ concentration are tabulated 

for each run in Table 7 . Our analyses of coupled chemistry and- 
diffusion in the post-vortex SST plume begins with the data tabu- 
lated in Table 6. 
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TABLE 5 a 


Plume centerline concentrations, [i] Q , variance of [HgO] 
[N 0 X ] lateral distributions, Cy, and integrated amounts 
of i th chemical species at t = 9 x 10 3 sec (2.5 hrs). 


v ,2 = 0, and = 1.6 x 10 -3 sec~\ 




Case 


Item* 

l 

2 

3 

4 

2 

a 

y 

1,87 x 10 3 

1.86 x 10 3 

1.83 x io 3 

1.83 x 10 3 


[0 3 ] o 

[KgOlo 

[H0 X ] O 

[N0] q 

[ho 2 ] 0 


[H 2 0) I 

[NO x ]i 

[NO] j 

[NOglj 

[HHOjlj 


2.1 x icr c 

7.6 x 10 15 

12 

U x 10 Xt - 

4.6 x 10 11 
6.0 x 10 11 


8.8 x 10 23 

1.3 x 10 20 

4.6 x 10 19 

7.7 x 10 19 


9.2 x 10 

7.6 x 10 15 

2.2 x 10 12 

1.4 x 10 x ^ 

7.8 x 10 11 


8,8 x 10 23 

2.5 x 10 20 

1.3 x 10 20 

1.2 x 10 20 

2 '-4 v 1 n 13 


1.2 x 10' 
7.7 x 10 ] 
1.1 x 10 1 
1.0 x 10 1 
7.9 x 10 ] 


8.8 x io 23 

1.2 x 10 21 

l.i x io 21 
1.6 x io 20 


7.2 x 10 

7.7 x 10 15 

2.2 x 10 13 
2.1 X IO 13 

9.6 x IO 11 


8.8 x 10 23 

2.5 x 10 21 

2.3 x 10 21 
2.0 x 10 20 

1.5 x IO 18 


ti] 0 = Plume centerline concentration, molecules/cm 3 

[i] = Integrated amount across the plume in a column of 

unit cross-section, molecules. 

2 o 

a y « Variance of [H g 0] and [NO ], m . 

The environmental values of [ 0 o 1 = 5 x 10 12 molecules/cm 3 : 

all other [il = 0. 3 e 

6 
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TABLE 5b 


Same as Table except v' 2 = 0 m 2 /sec 2 and 
j-|_* = 8 x 10"5 sec - "* - . 




Case 


Item* 

1 

O O 

w> 

4 

2 

c 

y 

1.87 x io 3 

1,87 x 10 3 

1,83 x 10 3 

1.83 x 10 3 

[0 3 ] o 

2.3 x 10 12 

1.3 x 10 12 

3..1 x 10 11 

11 

2.5 x 10 11 


[H 2 0] o 

^o 

[N03 o 

[N0 2 ] o 


[h 2 °]i 

[NO x ]i 


[HN0 5 ] I 


7.6 x 10 15 

l 2 

1.1 x io~ d 

7.6 x 10 11 

2.2 x 10 11 

8.8 x 10 23 


on 

1.2 X 10^ -' 


1 Q 

8.3 x 10~ J 

2.9 x 10 19 

1.4 x 10 19 


7.6 x 10 1J 

l 2 

2.2 x 10 

1.8 x 10 12 
3„Q x 10 11 

8., 8 x 10 23 

2.5 x 10 20 

1.9 x 10 20 

4.7 x 10 19 


1.2 x 10- 


7., 7 x 10 XJ 

1.1 x 10 13 

1.0 x 10 13 

1 1 

7.7 x 10"' 

8.8 x 10 23 

pi 

1 » 2 x 10 

1.2 x 10 21 


3 .8 x 10 J 


7.7 x 10^ 

2.2 x 10 13 

2.1 x 10 13 

6.8 x 10 12 

■ | 

8.8 x 10 23 
2.5 x 10 21 
2.4 x 10 21 


3.1 x 10 


[1] = Plume centerline concentration, molecules/cm 3 

[i]j = Integrated amount across the plume In a column of unit 
cross-section, molecules. 

a 2 = Variance of [H^O] and [N0 x ], m 2 . 

TP o 

The environmental values of [0,] = 5 x 10 molecules/cm- 5 ; all 

other [ i ] - 0. 8 e 
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TABLE 5c 



Same as 

Table 

5a 

— 2 

, except v 1 

= 0. 

01 

p 

m /sec 

2 

and 



V - 1 

.6 

X 10 

“5 

sec 

-1 















Case 






Item* 






2 



3 



!+ 


a 2 

y 

3.1 

X 

10 5 


3.1 

X 

10 8 

3.0 

X 

10 5 

3.0 

X 

10 5 

[0 3 ] o 

4.8 

X 

10 12 


4.7 

X 

10 12 

9.3 

X 

10 11 

1.5 

X 

10 12 

t H 2 °] 0 

6,0 

X 

10 14 


6.0 

X 

10 l4 

6.0 

X 

lO 1 ^ 

5.9 

X 

lO 1 ^ 

[ H0 x! o 

8,6 

X 

10 10 


1.7 

X 

10 11 

8.5 

X 

10 11 

1.7 

X 

10 12 

[N0] o 

2.1 

X 

10 10 


4,3 

X 

10 10 

5.3 

X 

10 11 

8.8 

X 

10 11 

[n° 2 ] 0 

6.2 

X 

10 10 


1.2 

X 

10 10 

3.2 

X 

10 U 

8.1 

X 

10 11 

H 

O 

OJ 

8.8 

X 

10 23 


8.8 

X 

10 23 

8.8 

X 

10 23 ' 

8.8 

X 

10 28 

[NQJ, 

1.3 

X 

<1 

10 20 


2.5 

X 

10 20 

1.3 

X 

10 21 

2.5 

X 

10 21 

[N0] I 

3.1 

X 

10^ 


6.4 

X 

10 19 

4.5 

X 

10 20 

1.1 

X 

10 21 

[HO^ 

9.2 

X 

lO 1 ? 


1.9 

X 

10 20 

8.0 

X 

10 20 

1.4 

X 

10 21 

t hno 3 3 -j- 

4.1 

X 

10 18 


4.2 

X 

0 

M 

00 

3.1 

X 

10 18 

2.8 

X 

10 18 


= Plume centerline concentration, molecules/cm 3 

[i]j = Integrated amount across the plume in a column of 
unit cross-section, molecules. 

a y = Variance of [HgO] and [NOJ, m 2 . 

The environmental values of [0 o ] . 5 x lo 12 molecules/cm 3 * 

all other [i] = 0 . 3 J e. J moiecuies/cm , 




TABLE 5d 

2 2 2 

Same as Table 5a, except v» = 0.01 m /sec and 

J x * =8 k 10"5 sec -1 . 



3.4 x 10 5 3.2 x 10^ 3.0 x lO'' 2.7 x 10 


4.7 x 10 J 

6.0 x 10 ] 

8.6 x 10 - 

4.4 x 10 J 

2.6 x 10 H 


4.6 x 10 J 

6.0 x 10 ] 

1.7 x 10 ] 
9.9 x 10 ] 

5.7 x 10 ] 


2.8 x 10 J 

6.0 x 10 3 

8.5 x 10 ] 

6.2 x 10 ] 


2.2 x 10 ' 


1.3 x 10 J 

6.4 x 10 ] 

1.8 x 10 ] 

1.6 x 10 3 

2.5 x 10 3 


8.9 x 10 23 8.9 x 10 23 8.8 x 10 23 8,8 x 10 23 


1.3 x 10 £ 


4.0 x 10 19 

2.1 x 10 19 


2.5 x 10 £ 


1.3 x 10 ‘ 


2.0 x 10 


1.4 x 10 


2.5 x 10 £ 


6.8 x 10 19 1.5 x 10 20 8.8 x 10 20 2.0 x 10‘ 

4.0 x 10 19 8.6 x 10 19 3.7 x 10 20 5.4 x 10‘ 


8.7 x 10 


[i] = Plume centerline concentration, molecules/cm 

[i], = Integrated amount across the plume in a column of 
unit cross-section, molecules. 

c 2 = Variance of [H^O] and [N0 X ], m 2 . 

The environmental values of [0 Q ] - 5 X 10 12 molecules/cm 

all other [± 1 _ = 0 , J 




TABLE 5e 


Same 

1 # = 
J 1 

as Table 5a, 

1.6 X 10" [i] * * * 5 

except 

sec”- 1 - 

v< 2 = 

0.02 m 2 /sec 2 and 






Case 



Item* 

1 


2 


3 

4 

a 2 

6.1 x 

10 5 

5.9 x 

10 9 



y 







[“A 

4.8 x 

10 12 

4.8 x 

10 1S 

N 

N 

O 

* — i 

0 

0J 

ixj 

1 » 

4.3 x 

10 14 

4,3 x 

10 14 

o 

0 

t NO x!o 

6.2 x 

10 10 

1.2 x 

10 U 

T 

T 

[M0] o 

1.5 x 

10 10 

3.0 x 

10 10 



[mo 2 ] 0 

4,5 x 

10 10 

9.0 x 

10 10 




8.8 x 

io 28 

8,9 x 

io 2 ' 8 

R 

R 

[NOJj 

1.3 x 

10 20 

2,5 x 

10 20 

U 

U 

[NO]! 

3.0 x 

io 19 

6.2 x 

10 19 

N 

N 

[NO^ 

9.2 x 

10 19 

1.9 x 

10 2G 



[hno 3 ] x 

4.2 x 

10 18 

4.2 -x 

10 18 




[i] = Plume centerline concentration, molecules/cm 8 

[ i ] -j- = Integrated amount across the plume in a column of 

unit cross-section, molecules. 

a 2 = Variance of [H o 0] and [NO 1 , m 2 . 

The environmental values of [0,] =* 5 X 10 12 moleeules/cm 8 : 

all other [i] - 0. 5 e 
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TABLE 5f 


2 

Same as Table 5a, except v 1 = 0.05 m^/sec 
j n * = 1.6 x 10“5 sec - - 1 -. 


1.5 x 10' 


1.4 x 10 


4. 

,9 

X 

10 A 

4. 

9 

X 

10' 

2 4 

,8 

X 

10 14 

2. 

8 

X 

10' 

4. 

.0 

X 

10 10 

7, 

9 

X 

10 

9. 

.3 

X 

o\ 

0 

i — ! 

1 . 

9 

X 

10 

2 , 

,9 

X 

0 

0 

5- 

8 

X 

10 

8, 

.9 

X 

10 2 3 

8, 

9 

X 

io : 

1, 

.3 

X 

10 20 

2. 

3 

X 

10 

3, 

,0 

X 

10^ 

6, 

2 

X 

10 

9, 

,2 

X 

10^ 

1, 

.9 

X 

io : 

4, 

.3 
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[1J = Plume centerline concentration, molecules/cm^ 

[ilj = Integrated amount across the plume in a column of 
unit cross-section, molecules. 

c 2 = Variance of [H Q 0] and [NO ], m 2 . 

j d X 

1 p 

The environmental values of [0_] = 5 x 10 molecules/cm 

all other [il = 0. ^ 







TABLE 6 


The effective eddy diffusivity, associated 

with turbulence levels chosen for post-vortex 
plume chemistry and diffusion calculations, based 
on = 2kt at t = 2.5 hrs. 

y 


Turbulence level 

Effective Eddy 

v’ (m /sec ) 

Diffusivity, ky 

O 

(m /sec) 

0 

0 

0.01 

1.7 x 10 1 

0.02 

3.2 x 10 1 

0.05 

8,2 x 10 1 





The fractional deficit of 0^ in a plume column, a, and at the pl ume centerline, 
b, as a function of the columnar N0 X , [NO x ]j, the turbulence v 1 and the 
photolysis rate, 

[NO, ] T Molecules 

X _L 

0^ 1.26xl0 20 2.52X10 20 1.26xl0 21 2.52x10 21 

v' 2 Def. = 

m 2 /sec 2 1.6xlO~ 5 8xl0~^ 1.6xlO~ 5 6xlO~ 5 1.6xlO~ 5 8xl0~ 5 1.6xlO” 5 8xlO -5 





The Fate of O3 in the Post -Vortex Plume 

The fate of 0^ during the first few hours of the SST plume 
life is of Interest for two reasons; l) the concentration of O3 

in the plume .determines the production of 0( 1 D) in the sunlit sky, 
and, therefore, the production of OH and, subsequently, HNO3; 

2) if the O3 is seriously depleted in the plume, the catalytic 
destruction of Oo by the N0 X cycle is much reduced (ref. 4). 
Depending upon how long this latter situation persists, the 0^ 
depletion may alter the estimate of the efficiency of the N0 X 
catalytic cycle in the general problem of environmental impact. 

In the first problem, i.e., the O3 role in HNO3 production, 
we must also be concerned about the depletion of 0^ in relation to 
the profiles of H 2 0 and N0 2 . Since [H 2 0] is a maximum at the plume 
centerline, and [N0 2 ] tends to be concentrated in the plume core 
(where [N0 X ] is a maximum), we have examined the plume centerline 
depletion of O3, as well as the columnar depletion across the plume. 

An inspection of Table 7 points up quickly under what circum- 
stances the depletion of 0^ becomes significant. The maximum 
depletion of O3 associates, not surprisingly, with minimum turbu- 
lent mixing (O3 is not entrained from the environment ) , and with 
maximum values of NO . However, in the sunlit sky, the regenera- 
tive capacity for the production of O3 in the plume is sufficient 

to limit the columnar depletion of O3 to about 40 to 50 per cent, 
although the centerline concentration of O3 has essentially gone to 
zero when [N0 x ]j exceeds 3 x 10 21 molecules. (The importance of 
this O3 depletion to HNO3 formation is discussed in the next 
section.) The general depletion of O3 in the plume appears to 
become significant when the turbulence level is less than 0.01m 2 /sec 
and the [NOjj exceeds 1.26 x 10 21 molecules. The latter number 
associates with an emission index for' N0 X of about 0.005, which is 
only about l/40th current N0 X emission control technology (ref. 19). 
Since near-zero levels of small -scale turbulence may also be more 
the rule than the exception in the lower stratosphere, we conclude 
that the early depletion of 0^ in the SST plume cannot be dismiosed 

out of hand. 

Within the ranges of the input variables tested, the data in 
Table 7 permit an "engineering" estimate of the rela tion ships 

between 0 3 depletion and the controlling variables, v 1 and [N0 x ]j. 
First, a plot of plume centerline O3 deficit, AO^Gjj ), vs the 

columnar deficit, A0 3 (Col), (Fig. 18) shows these quantities are 

highly correlated for all values of v ,£1 and [H0 x ]-j- tested. This 
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AO3 (Col), Fractional O3 deficit 
in plume column 



AO 3 (q_ ), Fractional 0 3 deficit 
on plume (J_ 


Figure 18. The relationship of 0^ deficits at the plume 

centerline to the columnar CU deficit at 
t = 2.5 hrs . 5 


result suggests that, for this wide variety of conditions, the O 3 
concentration profiles have a high degree of self -similarity . 
Within the ranges tested. 


a0 3 ( c l ) 

ao^TcoT) 


= 4 


( 66 ) 


The partial dependence of AO^Cj^ ) on [ N0 X ] j and v 1 is shown 
in Fig. 19, These data yield the relationship, 


[N° x ] 

AO (^ ) = C -=4 
5 1 ( 7 * 2)2 


(v ' 2 > . 001 ) 

([N0 X ] I < 2.5 x 10 21 ) 


(67) 


where C-^ - 3 x lO" 2 ^. Since the chemistry calculations have shown 
that [0 3 ) is insensitive to [H 2 0] until [H 2 0] exceeds 10 17 mole- 

cules/cm3 (an excessive value in the post-vortex plume), Eq. (67) 
captures the essential dependence of the plume 0 ^ deficit on plume 
and turbulence parameters. We note in particular that a signifi- 
cant 0 ^-hole develops at the core of the plume when the ratio 
3 'o L pp 

[N° x ]i/(v' ) 2 exceeds 3 x 10 , a value only slightly beyond the 

range we have tested. 


The Fate of NO in the Post-Vortex Plume 
.X 

in examining the predicted fate of N 0 X in the early post- 
vortex plume, we are primarily interested in the rapid formation, 
of HNOo. This interest, of course, resides in the strong buffering 
of the NO catalytic cycle for the destruction of which conver- 
sion to HNO^ represents. The chemistry calculations have shown 
this conversion is rapid only when high concentrations of H 2 0 and 
Oo are associated with significant concentrations of N0 2 . If this 
conversion of N0 2 to HNO 3 is to be a significant factor in the SST 
problem, it must occur before the H 2 0 is diluted by turbulent 
mixing. These early plume calculations are, therefore, particularly 
pertinent to this phase of the problem. 
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("b) £ ov 



[NOjj molecules 


Figure 19. The partial dependence of the 0^ centerline 

deficit at 2.5 hrs on 1) the level of turbu- 
2 

lence, v* , and 2) the amount of NO in the 
plume, [N0 X ] I . X 
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The results of the model calculations of HNO 3 formation during 

2.5 hrs and for the various conditions listed in Tables 3 and 4 
are tabulated in Tables 8 and 9 and displayed graphically in Fig. 20. 
Inspection of either of these displays reveals very quickly that the 
fraction of N0 X converted to HNO 3 is greater than 10 per cent only 

for minimum values of [N0 x ]j and maximum values of jj* . If the 


conditions tested are at all representative of the SST post-vortex 
plume, the early formation of HNO 3 cannot be reckoned as a signifi- 
cant process. This conclusion would be modified foi cases where 
j * and [HpO] are somewhat higher than the limits used here, i.e., 

higher in the stratosphere and for engines with larger fuel consump- 
tion rates. But, for the present problem, it appears early HNO 3 
formation provides only a minor safeguard against the later 
catalytic destruction of O 3 by NO and NOg* 

It is interesting to note that the formation ol HNO 3 is 
enhanced by turbulent mixing. The reason for this behavior is, of 
course, ascribable to the fact that the entrainment of environmental 
0 o more than offsets the dil ution of the other reactants. This 

general Interplay among v' 2 , O 3 , and N0 X In the formation of HNO 3 
can be seen when we plot the columnar O 3 deficit against [HNO^Jj/ 

[NO ] T , as shown in Fig. 21. This figure only identifies the 
inverse association of the 0. deficit with the suppression of HNO 3 
formation as a function of the turbulent mixing. However, this 
compar 1 son^oint s up the mutual fates of N 0 X and O 3 namely, when 
0o is depleted, due mainly to increasing N0 X concentrations, the 
formation of HNOo is suppressed. Although HNO 3 formation does no 

now appear significant, we may note in passing that any reduction 
?n N0 P em?L?ons not only reduces the total N0 X in the stratosphere, 
but also leads to a larger fraction of the N0 X being converted 
early -on to HNOn. It remains to be seen as to whether this compou 
ing effect of N 0 X reduction produces a significant change m N0-N0 2 
concentrations, and therefore In O 3 destruction, at uselul opera- 


tional limits. 


CONCLUSIONS 


The final evaluation and interpretation of these simulations 
of chemistry and diffusion in highly perturbed stratospheres mu 
be deferred^ until validation data regarding the turbulence field 
and the chemical composition of the early SST ex haust Pi™ 
been obtained experimentally. This final prooi o y ■ 

"model becomes particularly important t because^! the 

unusual sensitivity of the xares 01 U 3 ailu ,v x A 
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TABLE 8 

The total formation of HNCu in 2.5 hrs for various Initial 

5 1 A 

conditions and turbulence levels, [HNCu] T in 10 molecules 



[NO ] T Molecules 

X 


Turbulence 

level 1 . 26 xl 0 2 

2.52xl0 20 1.26 x10 21 

2.52x10 

v ' 2 J x * = 

m 2 /sec 2 l.bxlO - ' 5 8xl0 - ^ 

1.6xlO -5 8xl0 -5 1.6xlO" 5 8xl0 -5 

1.6xlO -5 8x 


0 

2.87 

14.4 

2.33 

11.8 

1.08 

1.08 

1.49 

3-1 

.01 

4.13 

20.6 

4.i6 

20.3 

3.10 

i4.o 

• 2.73 

8.7 

.02 

4.20 


4.26 






.05 

4.30 


4.35 









The fraction of NO 

x 


converted to HNO 


3 


in 2.5 hrs . 


[ NO x ] I Molecules 


Turbulence 

level 1.26xl0 20 2.52xl0 20 1.26xl0 21 2.52xl0 21 


m 2 /sec 2 1.6x10 5 8x10 5 1.6xlO~ 5 8x10 5 1.6x10 5 8x10 5 1.6xl0 -5 8x10 5 

0 .023 .114 .009 .047 <.001 .001 <.001 .001 

.01 .035 .164 .016 .081 .002 .011 .001 .004 

.02 .033 .017 

.034 


05 . 


.017 





CO 

<D 

o 

O) 

o 
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Figure 20. The dependence of fractional conversion of 
N0 X to HNO^ during 2,5 hours on the amount 

of NO in the plume and two values of the 
photolysis constants. 



A0 3 (Col) 



.gure 21. Joint values of plume 0^ depletion 
and HNO^ formation in 2.5 hours as 

calculated “by 'the coupled diffusio: 
chemistry model. 




in the plume. Best estimates of this quantity now place the N0 X 
emissions and initial plume N0 X concentrations precisely in this 
range of greatest sensitivity, i.e., 10 -*- 2 to 10 ^ molecules/cm 3 
in the core of the post-vortex plume, and generally In the range 
10 10 to 10 12 in the other parts of the plume. The remarkable 
decline of 0, concentrations as the N 0 X concentrations go above 
1012 molecules/cm 3 poses a "continental divide" effect insofar as 
the fates of 0^ and HNO^ formation in the plume are concerned. When 
NO concentrations go above 10 ^ 2 molecules/cm 3 , O 3 concentrations 
are rapidly suppressed, and HNO^ formation is inhibited due to the 
lack of 0( 1 D) . When [N0 X ] exceeds 10 13 the destruction of O 3 is 
essentially complete , resulting in an "O^-hole" in that region 
(generally the core) of the plume. The importance of this 0^ 
deficit to other chemical reactions, such as the fate of unburned 
hydrocarbons, has not been investigated in the present study. 

When [NO ] < 10 12 molecules/cm 3 , O 3 is only very slowly depleted 

and, given H 2 0 concentrations above lO 1 ^ to 10 13 molecules/cm 3 , the 
formation of HNO 3 is predicted. However, on the basis of plume 
conditions examined here, significant conversion of N0 X to HNO 3 
in a time period of a few hours occurs only when [HgO] is present 
in excess of about lO ^"^ 1 mo lecules/ cm 3 and [N0 X ] is less than 10 
molecules/cm 3 . These are incompatible requirements for present 
H 2 0 and N0 X emission estimates, since they require a minimum dilu- 
tion of N0 X of about two to three orders of magnitude, but a 
maximum dilution of H 2 0 of only ten. Since H 2 0 and N0 X emanate 
from a common source, differential dilution of these exhaust, 
constituents in the ratio of one to two orders of magnitude. is 
not physically possible. On this basis, we conclude that.signifi 
cant HNOq formation in the first few hours of plume life is not 
likely under the emission estimates presently used for H 2 0 and N0 X . 


We do note, however, that these simulations provide a basis for 
estimating the change in the environmental impact of reducing SST 
NO emissions which go beyond the simple proportional reduction in 
N 0 X and N0 2 . If the initial post-vortex concentrations of N0 X are 
reduced below 10 12 molecules/cm 3 ^by emission controls, but the H 2 0 

concentrations remain in the 10 l6 molecules/cm 3 range, then in the 
sunlit sky, a significant fraction of the N0 X can be converted to 
HNOo within a few hours. Since the photolysis oi HNO 3 , to recover 
NO ~*and OH, is slow, this conversion, if realized, would effective- 
lv 2 buffer the N0-N0 2 catalytic cycle for the destruction of 0^. 

Tust where these compounding effects of proportional reduction of 
[NO] and [NOg] by emissions controls and HNO^ formation sets m 
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depends upon reliable estimates of N0 X emission rates and turbulent 

mixing in the post-vortex plume. But, the simulation models for 
evaluating the payoff of N 0 X emission controls in decreased destruc- 
tion of 0 ^ are available. 

Finally, the role of stratospheric turbulence in the determina- 
tion of the fates of 0^ and HNCb are now much clearer. The primary 

roles of turbulent mixing are 1 ) to entrain environmental 0 ^ into 
the post-vortex plume, and 2 ) to dilute the N 0 X and H 2 0 concentra- 
tions in the plume, and, of course, to spread these exhaust 
constituents through an increasing volume of the stratosphere. We 
have argued that near-zero values of small-scale atmospheric turbu- 
lence should be quite common, indeed, the rule in the lower strato- 
sphere. Under these circumstances the mixing of exhaust materials 
with the environmental 0 -^, and the dilution of exhaust constituents 

are both very slow, and the chemical fate of exhaust materials is 
controlled by their relative abundance within the plume. During 
the time this nonturbulent condition persists, plume segments remain 
more or less intact (they may be distorted but not diluted by larger 
scale wind shear), and the internal balances of [0^], [NO], [N0 2 ], 

and [HNO^] operate on a "do-nothing" cycle. This very real possi- 
bility becomes significant to the environmental impact assessment 
of SST operations when these intact plume segments persist for a 
significant fraction of the residence time of exhaust materials in 
the lower stratosphere . 

At this point the lack of hard data on the turbulence history 
of individual small volumes of air in the lower stratosphere 
precludes a definitive assessment of this effect. What is needed 
are measurements of the time history of turbulence in the vicinity 
of free-floating balloons which track these individual volumes of 
air as they circumnavigate the globe. With appropriate turbulence 
instrumentation (the primary requirements are for resolution of 
vertical and lateral gusts down to rms values of the order 2-5 cm/ 
sec), the basic problems of turbulent mixing, the sporadic nature 
of this mixing, and the recurrence period for significant turbu- 
lence could be resolved. These data would permit much more realistic 
and reliable estimates of the fates of SST exhaust materials and of 
environmental 0 -^ than are now available. 

The problem of the residence time of SST exhaust materials in 
the lower stratosphere, given the highly specialized a ltitude and 
geographical locations oT the 80 T flight routes , goes beyond the 
immediate purposes" of the present study somewhat, but is repeated 
here for emphasis in interpreting the significance of early plume 
chemistry and diffusion. The validity of the assumption of a 
laterally well-mixed lower stratosphere, in which the contribution 
of SST exhausts are treated as a fractional perturbation of the 
background values of N0 X and H o 0, hinges on a reliable estimate of 
this residence time. eL 
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APPENDIX A 


The Effect of Incomplete Mixing During Vortex 
Decay on the Conversion of NO to NO^ 

We examine the effect of neglecting the mixedness term in the 
vortex decay calculations. Since the primary reaction of interest 
during the vortex decay period is the conversion of NO to NO 2 by 
O 3 , we can examine this effect by using an earlier version of the 
coupled model in which HgO was neglected, but the term [O 3 ] 1 [NO ] ’ 

was retained in the calculation ( ref. 5). In this case the 
mean time rate of change of [NO] is given by 

dlNOl = _ d_ v ,^ NO y, , k 2 ([N0][0 3 ] + [NO] • [0 3 ] . ) + j 2 [N0 2 ] (Al) 


Then 

00 

/M01 dy 
- 00 


J V [NO] ' dy - k 2 j ( [NO] [Og] + [ NO ] • [ 0 3 ] 1 ) dy 


r 

/ f r \T„ ' -i - TT’ 


00 

r\ 

+ Jo J [NO 2^ dy ( A 2 ) 

-00 


But 


00 


r d[N0~ 

j dt 


dy 


00 

n 


= ft / [ »] dy 


d [NO ] j- 
dt 


(A3) 
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b 


' [ NO ] 1 dy = 0 


(A4) 


and 


3 2 j [N0 2 ] dy = J 2 [N0 2 ] 1 

—bo 

Then 


(A5) 


d[NO] 

dt 


T n _ 

-= - k 2 / ( [ NO ] [ 0 ] + [ NO ] * [ 0 ]«) dy 


+ J 2 [ NO ^ ] i 


If the mixedness term [ NO ] 1 [ O^ ] ’ is set equal to zero, 


(A6) 


dlNO)/ 

dt 


r _ 

k 2 j [NO] [0^] dy + j 2 [N0 2 ] J 


(AT) 


and 


d[NO] x /dt - j 2 [TJ0] I [ [NO] [0 3 ] + [N0]‘[0 3 ]-] I 


d[NO] x /dt - j Q [ NO ] 


[ [NO] [0 3 ] ] x 


(A8) 


To estimate the magnitude of the term on the right-han d s ide 

2 

of Eq. (A8), this model was run for the vortex decay case, v' 

0.1 m 2 /sec2, [O^Iq = 4 x lO"*"^ molecules/cm^, [N0] o = 10^3 mole- 
cule s/cm3, [N0 2 ] = 0, and c2 q = 100 m 2 . The values of 
[ [NO] ' [ 0 3 ] 1 ] j/( [NO] [ 0 3 ] ) j for 400 seconds are shown in Pig. A-l. 
As can be seen there, the effect of mixedness is to suppress the 
reaction between [0^] and [NO], This effect reaches a maximum of 

about 30 per cent around 75 seconds and has decayed to less than 
10 per cent after 300 seconds. On this basis the coupled model 
used in the present study overestimates the conversion of NO to 
NOg by approximately 20 per cent during the 200 seconds assumed 

for vortex decay turbulence. 
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Time, seconds 


Figure Al. Time history of the mixedness, integrated 
over the plume width, as a function of 
turbulence and stability. 
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